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29 OCTOBER 1945 


IN DIFFERENTIAL ANALYZER ROOM. 


H. B. V. Busn, H. L. Hazen, S. H. 


Conference on Advanced Computation 
Techniques 


During the recent war there was tremendous development of certain 
types of computing devices. The big Harvard-IBM machine indicates the 
way in which a whole battery of IBM equipment can be brought under 
automatic electric control. Electronic devices have been constructed, or 
are under development, which promise astronomical speeds for numerical 
processes. At the Bell Laboratories there has been important development 
of automatic relay computing. Thus, with war’s end these and other similar 
developments suggest that there will soon be available mechanical and 
electrical computing equipment which, in terms of speed and flexibility, 
will completely outdistance anything thought of before. 

With this in mind the desirability of a Conference of those who, in this 
country and in England, have been recently directly and indirectly con- 
nected with the output of large computing machines became evident, and 
the National Research Council Committee on Mathematical Tables and 
Other Aids to Computation took steps for its organization. Most of the 
responsibility of this organization was placed in the hands of the Com- 
mittee’s Subcommittee Z (On Calculating Machines and Mechanical Com- 
putation) which then! had the following membership: 


L. J. Comrir, Chairman J. C. P. MILLer 
S. H. CALDWELL, Vice-Chairman G. R. Stisitz 

H. H. AIKEN I. A. TRAVIS 

D. H. LEHMER J. R. WoMERSLEY. 


Since the chairman had not yet arrived in this country, the vice-chairman 
carried out most of the detailed arrangements. 

The Massachusetts Institute of Technology not only offered its facilities 
for such a gathering but also suggested that it should coincide with the first 
public demonstration of its new Differential Analyzer. The dates for the 
gathering were fixed as October 29, 30, 31, 1945. There were the following 84 
guests who accepted our invitation: ° 


Name Representing 

Abramowitz, Milton............. Math. Tables Project, 150 Nassau St., New York City 

Aiken, H. H. (Commander)........ Computation Project, Cruft Laboratory, Harvard Uni- 
versity 

Ground, Maryland 

Archibald, R. C. (Prof.).......... Nat. Res. Council Comm. (Brown University) 

Arnold, H. A. (Lt. Commander). .. .Computation Project, Cruft Laboratory 

Bennot, Maude (Miss)........... 25 Ridgeway Avenue, Needham, Mass. 

Bergman, Stefan (Dr.)........... Brown University 

ok Computation Project, Cruft Laboratory 

Computation Project, Cruft Laboratory 

Bell Telephone Laboratories, New York City 

Brainerd, J. G. (Prof.)........... University of Pennsylvania 

Brendle, Ruth (Miss)............ Computation Project, Cruft Laboratory 
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Name Representing 

EE Office of Research & Inventions, 1 Park Avenue, New 
York City 

Bush, Vannevar (Dr.)............ Carnegie Institution, Washington, D. C. 

Capewen, S. Bi. (Prof.)........... Nat. Res. Council Comm. (Mass. Inst. of Technology) 

Campbell, R. V. D. (Ensign)....... Computation Project, Cruft Laboratory 

General Electric Co., Schenectady, N. Y. 

ae General Electric Co. 

OE) Nat. Res. Council Comm. (London, England) 

Concordia, Charles.............. General Electric Co. 

Terre Office of Research & Inventions 

Cunningham, L. E. (Dr.)......... Ballistic Research Laboratory 

SS) Ballistic Research Laboratory 

6:5 National Cash Régister Co., Dayton, Ohio 

International Business Machines Corp., 590 Madison 
Ave., New York City 

Sere University of Pennsylvania 

re Mass. Inst. of Technology 

Goldstine, H. H. (Capt.).......... Ballistic Research Laboratory 

SS ee Office of Research & Inventions 

International Business Machines Corp. 

Haskins, Robert (Lt.)............ Office of Research & Inventions 

Mass. Inst. of Technology 

rrr Mass. Inst. of Technology 

Hopper, Grace (Lt.)............. Computation Project, Cruft Laboratory 

Mass. Inst. of Technology 

8 Institute Math. Statistics (Harvard University) 

(Dr.)............. Mass. Inst. of Technology 

rere General Electric Company 

Nat. Res. Council Comm. (Ballistic Research Labora- 
tory) 

EE ee Mass. Inst. of Technology 

Math. Tables Project 

Eastman Kodak Co., Rochester, N. Y. 

International Machines Corp. 

General Electric Co. 

Mauchly, J. W. (Prof.)........... University of Pennsylvania 

Mayall, Margaret (Mrs.)......... Harvard Observatory, Harvard University 

Office of Research & Inventions 

National Cash Register Co. 

Myers, J. E. (Ensign)............ Naval Research Laboratory, Anacostia, Md. 

von Neumann, John (Prof.)........ Institute Math. Statistics (Institute for Advanced 
Study) 

Nyquist, Harry (Dr.)............ Bell Telephone Laboratories 

re Foxboro Instrument Company, Foxboro, Mass. 

8 rere Applied Mathematics Panel (350 Fifth Avenue, New 
York City) 


C. Office of Research & Inventions 
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Name Representing 
Sandomire, Marion M. (Mrs.)..... Bureau of Ships, Washington, D. C. 
Naval Research Laboratory 
Bell Telephone Laboratories 
Sternberg, Sidney (Lt.)........... Office of Research & Inventions 
8 eee University of Vermont, Burlington, Vt. 
Taylor, Richard (Prof.)........... Mass. Inst. of Technology 
Thomas, G. B. jr. (Prof.)......... Mass. Inst. of Technology 
errr Eastman Kodak Co., Rochester, N. Y. 
Mass. Inst. of Technology 
We Office of Research & Inventions 
Wallman, Henry (Dr.)........... Mass. Inst. of Technology 
Warren, S. R. ir. (Paol.). .......4.. University of Pennsylvania 
Watson, F. G. (Lt. Commander). ... Hydrographic Office, Navy Dept., Washington, D. C. 
Wickinad, H. P. (it.)............ Office of Research & Inventions 
Wiener, Norbert (Prof.).......... Mass. Inst. of Technology 
re National Cash Register Co. 
tS Chance-Vought Aircraft Div. of the United Aircraft 

Corp., Stamford, Conn. 

Wundheiler, A. W. (Dr.).......... Bureau of Ordnance 


Most of these delegates arrived on October 29 and attended demonstra- 
tions of the working of the new Differential Analyzer, of which there is a 
detailed illustrated description by VANNEVAR Buss and S. H. CALDWELL,” 
in Franklin Institute, J., v. 240, Oct. 1945, p. 255-326. Our frontispiece 
reproduces a photograph taken at the opening session of the Differential 
Analyzer demonstration. Grouped about the Main Control Panei are (3) 
Mr. Hazen, head of the Department of Electrical Engineering at M.I.T., 
(2 & 4) Messrs. Bush & Caldwell, who were initially responsible for develop- 
ing the Analyzer, and (1) Mr. Phillips, head of the Department of Mathe- 
matics at M.I.T. In the evening the Institute invited the Delegates to a 
delightful dinner at the Engineers Club on Commonwealth Ave. The 
addresses of JAMES R. KILLIAN, JR., executive vice-president of the Insti- 
tute, presiding officer, GEORGE R. HARRISON, dean of science at the Institute, 
LEsLIE J. Comrig, chairman of Section Z, and VANNEVAR Busi, director 
of the Carnegie Institution of Washington, were uniformly of exceptionally 
happy content. 

The Program carried out on Tuesday, Wednesday, October 30, 31, with 
L. J. Comrie and S. H. CALDWELL presiding, was as follows: 


Tuesday Morning, M.I.T. 
Address of Welcome to the Institute by President K. T. Compton. 
Introductory Remarks by VANNEVAR BusH. 
“‘Some Industrial Applications of Machine Computing Methods” by F. J. MAGINNIss. 
“Numerical Computation by the Use of Telephone Relays” by S. B. WILLIAMs. 
Tuesday Afternoon, M.1.T. 
“Some Principles of Electronic High Speed Computing” by J. G. BratnerD and 
J. P. Eckert, JR. 
“A General Summary and a Look at the Future” by J. von NEUMANN. 
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Tuesday Evening, Cruft Laboratory, Harvard University 

“The Automatic Sequence Calculator’’® by H. H. AIKEN, and members of his staff. 
Wednesday Morning, M.I.T. 

“Simulative Techniques as an Aid to Analysis” by G. A. PHILBRICK. 

“Application of Punched-Card Methods to Scientific Computation” by W. J. Eckert, 

presented by J. C. McPHERSON. 

“The M.I.T. Center-of-Analysis Research Program” by S. H. CALDWELL. 
Wednesday Afternoon, M.1.T. 

“The Scientific Application of Commercial Calculating Machines” by L. J. Comrig. 


The Conference was most notably successful, and one heard on every 
side expressions of the hope that such a Conference might become an annual 
event. 

The Committee on Mathematical Tables and Other Aids to Computa- 
tion is indeed greatly indebted to officers of the administration and members 
of the staff of the Massachusetts Institute of Technology who did every- 
thing in their power to provide a setting for securing fruitful results from 
such an assembly of experts. 

R. C. A. 
1 Beginning with Nov. 1, 1945, S. H. C. became chairman, and L. J. C. vice-chairman. 


2 See MAC 19. 
* For illustrated notes on this Calculator see MAC 20. 
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275{A, E, F]—PEDER PEDERSON,, (a) Uber die numerische Berechnung der 
Kettenbriiche nebst einer Berechnung der Grundzahl der natiirlichen Loga- 
rithmen, 36 p., 1940; (b) Berechnung der Grundzahl e der natiirlichen 
Logarithmen mii 606 Dezimalen, 17 p., 1942; (c) Fortsetzung der Berech- 
nung der Grundzahl e der natiirlichen Logarithmen bis zur 808. Dezimal- 
stelle, 21 p., 1944. Denmark, Geodetisk Institut, Meddelelse, nos. 14, 
16, and 17 respectively. 14.322 cm. 


These three papers are accounts of the author’s calculations of the fundamental con- 
stant e. Values of e are given to 404D in (a), to 606D in (b) and to 808D in (c). The method 
of calculation was based on the continued fraction 


2 1+6+10+14+--- 


which is due to EuLEr,! but is applicable to any regular continued fraction. The general 
idea is set forth in (a) and may be roughly described as follows. Let x be the exact value 
of a continued fraction and let x, be a given approximation to x which is correct to m deci- 
mals. Let also A;/B; be any convergent of the continued fraction in which k is large enough 
so that the value of this convergent would give x to more than m decimals. This means 
that, roughly speaking, B, has about n/2 digits. Then x,B; will be a decimal whose first 
n/2 significant figures are the digits of A; and whose next /2 digits are zeros. Therefore 
if we compute denominators B; only, and compute only the digits of the product xB: 
which lie near the mth significant figure, the appearance of zeros (or nines) at this place 
serves as a check on the digits of x,. Moreover, the method may be used to extend the 
value of x beyond the mth decimal place by simply filling in those digits which will produce 
more zeros until the accuracy implicit in a particular B; is exceeded. 

The author started with SHanxs’ 137D value of e verified by GLAISHER,? and went 
on to check the value of Shanks* (205D, last 18 wrong) and of Boorman‘ (346D, last 
123 wrong). The 223 correct decimals of Boorman was then the basis of an extension to 
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404D, and then to 606D. At this point the author became aware of a 707D value given in 
1926 by the reviewer,’ which he verified and used to extend the value of ¢ to the impressive 
result of 808D. The suggested method of checking and extending the value of a continued 
fraction, given in 1926, is thus independently discovered and greatly elaborated by the 
author. 

The determination of the digits of e, one at a time, seems to suggest that the author 
used no calculating machine to perform this immense calculation. 

Besides the values of e mentioned above the author gives all the denominators of the 
convergents up to By, a 408 digit integer. These interesting numbers begin with B, = 1, 
B, = 7, B; = 71, By = 1001 and were found by the recurrence formula 

Bix = (4k + 2)Be + Bea. 
They have interesting divisibility properties discovered by D. N. Lenmer® and further 
studied by the reviewer.? The author gives no description of this part of the work and of 
the checks which are available to insure the accuracy of this, the most hazardous part 
of the calculation. 

Since the periodical in which this value of e to 808D appears is highly inaccessible in 
this country, it seemed advisable to reproduce it below. 


e = 2.71828 18284 59045 23536 02874 71352 66249 77572 47093 69995 
95749 66967 62772 40766 30353 54759 45713 82178 52516 64274 
27466 39193 20030 59921 81741 35966 29043 57290 03342 95260 
59563 07381 32328 62794 34907 63233 82988 07531 95251 01901 
15738 34187 93070 21540 89149 93488 41675 09244 76146 06680 
82264 80016 84774 11853 74234 54424 37107 53907 77449 92069 
55170 27618 38606 26133 13845 83000 75204 49338 26560 29760 
67371 13200 70932 87091 27443 74704 72306 96977 20931 01416 
92836 81902 55151 08657 46377 21112 52389 78442 50569 53696 
77078 54499 69967 94686 44549 05987 93163 68892 30098 79312 
77361 78215 42499 92295 76351 48220 82698 95193 66803 31825 
28869 39849 64651 05820 93923 98294 88793 32036 25094 43117 
30123 81970 68416 14039 70198 37679 32068 32823 76464 80429 
53118 02328 78250 98194 55815 30175 67173 61332 06981 12509 
96181 88159 30416 90351 59888 85193 45807 27386 67385 89422 
87922 84998 92086 80582 57492 79610 48419 84443 63463 24496 


D. H. L. 


1L. Euter, Introductio in Analysin Infinitorum, v. 1, Lausanne, 1748, p. 319; Opera 
Omnia, 8, Leipzig and Berlin, 1922, p. 388. 

28, Oy A L. GLAISHER, B.A.AS., Report, i871, p. 16-18 (Sect. trans.), also Report, 1873, 
p. 1 i 


2W. SHANKS, Royal So. London, Proc., v. 6, 1854, p. 397. 

‘J. M. Boorman, Math. Mag., v. 1, 1884, p. 204. "another calculation was made by 
Fr. Tichaneck in 1893 to 225D (correct to within a unit in the last place), Jahrbuch i. Ff 
Fortschritte d. Mathem., v. 25, p. 736; there is here a misprint in the 43rd decimal place, 
for 6 read 0. 

5D. H. Lenmer, Amer. J. Math., v. 48, 1926, p. 139-143. 

6D. N. Lenmer, Amer. J. Math., v. 40, 1918, p. 375-390. Also Nat. Acad. Sci., Proc., 
v. 4, 1918, p. 214-218, 

7D. H. Lenmer, Annals Math., v. 23, 1932, p. 143-150. 


EpiroriAL Note: For various values of e*, see MTAC, v. 1, p. 54-55. 


276[B, K, M].—A. N. Kotmocorov, (a) “Chislo popadanii pri neskol’kikh 
vystrelakh i obshchie prinfSipy ofSenki effektivnosti sistemy strel’by” 
(The number of hits in several rounds, and the general principles of 
estimating effectiveness of a system of firing); (b) ‘“Iskusstvennoe 
rasseivanie v sluchae porazhenita odnim popadaniem i rasseivaniia v 
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odnom izmerenii” (Intentional dispersions in case one hit is fatal anu 
dispersions in one dimension), Akademifa Nauk SSSR, Leningrad, 
Matematicheskii Institut imeni V. A. Steklova, Trudy, v. 12, June 194.- 
p. 7-45, 95-106. 


This issue of the Trudy is devoted to a collection of articles on the theory of gunnery, 


edited by Kolmogorov. There are four tables (p. 96-106). Table 1 gives values of 
a? 
H(m,a)=1-—€e (145454 oe 


for a = [.1(.1)15; 5D], m = 1(1)11. 
For the same range T. 2 gives the values of 


V?H(m, a) = H(m, a) — 2H(m — 1,a) + H(m — 2, a). 


In T. 3, values are given for : 
r(u) = = fin {1 + u) — 


for u = [0(.05)1(1)5; 3D], with A; (z) = = i edu, p = 476936... 

T. 4 gives the values of P(N) = (so) — 2sop(so), o(z) = so(N) = Ni, 
and d?(N) = 450%, for N = [0(.1)1(.2)5(.5)10; 3S] with A. 

On p. 104-106 are graphs of (a) H(m,a), m = 1(1)10, a = 0(1)15; (b) R(k, c, n), 
k = 3, c = 4/15, m = 2(2)6, 12(6)24, 0<a<4.8, the values of R being given in a table 
on p. 33 for given values of m and a; (c) R, k = 3, c = 4/15, n = 3, 6(6)24, 0<a<2.4, 
the values of R being given on p. 34 for values of m and a. 


R. C. A. 


277(D, J).—Great Britain, H.M. NavuticaL ALMANAC OFFICE, Tables for 
the Summation of Trigonometrical Series. Department of Scientific Re- 
search and Experiment, Admiralty Computing Service, November, 1945. | 
No. SRE/ACS. 93. 24 leaves mimeographed on one side of the paper. 
20.2 X 32.9 cm. This publication is available only to certain Government 
agencies and activities. 


These tables give the values of sin mx and cos mx for nm = 1(1)10, x = 0(1°)180°, and 
also for radian arguments x = 0(.01)3.14. The values corresponding to m = 1, 2, 3, and 4 
are given to 5D; those corresponding to m = 5, 6, and 7 to 4D; and those corresponding to 
n = 8,9, and 10 to 3D. 

These tables were compiled for use in the summation of trigonometric series of the 
form 2(A, cos nx + B, sin nx) when the convergence is fast enough for the sum of the 
first ten terms to give adequate accuracy. In conjunction with a table of powers, the tables 
may also be used to convert polar into cartesian coordinates. 

ARNOLD N. Lowan 


EprroriaL Note: For radian argument the above-mentioned table ranges, at varying 
intervals, from 0 to 31.4 radians, 5-3D. The NYMTP Tables of Sines and Cosines for Radian 
paar yiy M00; SB” among other tables, sin x and cos x, for x = [0(.001)25; 8D], and 
1 q 


278[E, M].—Louis Dunoyer, “‘Calculs sur l’émission en bout des tubes 
lumineux et notamment des tubes électro- ou photoluminescents,” 
Revue d’Optique, v. 18, 1939, p. 235-241. 1624.2 cm. 

There are several tables and graphs. The tables are as follows: 


(a) y = e~*/8i9 7 sin y, fora = .5, and a/sin y = [.5(.01).53, .55(.5).75, .85, 1(.25)1.5(.5)3(1)6; 
5D}. 
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tb) H1(m) = f°" sin yar, for = 2°(4°)90°; SDI, and a = .5, m = cot u. 
\2) H(m) = sin dy, for a = 0, .2,.5, 1, 2, and m = (0, 01(.02).99, 1.05(.1)1.95, 

also .1(.2)8.5, 9.5; 4D]. 
(4) + H(M = f° sin ada, 

H(M — m) = sin cot u’ = M — m, for 

M = (0(.1)1(.2)4(1)10, ©; 4D] and a = k = 0, .2, .5, 1, 2. 

There are graphs of H(m), a = k = 0, 2, 5, 1,2,and0 < m < 5;0f Hm + H(M — m) 
for M = .2(.2)1(.5)2(1)5,0<m< 5; andof [H(m) + H(M —m)] A(m)dm for a =k = 0, 
.2, 5, 1, 2, and 0 < _m 4. 
R. C. A. & S. A. Jorre 


279[F].__N. G. W. H. BEEGER, “Extension of the table of least exponents £ 
for which 2§=1 (mod p).”” Nieuw Archief v. Wiskunde, s. 2, v. 20, 1940, 
p. 307-8. 15.7 X 23.6 cm. 


The number ¢ referred to in the title is a function of p familiar to number theorists 
and usually called the exponent to which 2 belongs (mod ). It is known that ¢ is some 
divisor of p — 1. The other factor » = (p — 1)/é has been called by A. J. C. CunnincHam 
the residue index of 2 (mod p) and this function, because it is nearly always much smaller 
than », is usually tabulated in lieu of ¢. The present note contains a one-page table giving » 
for each prime p between 300000 and 310000. As a matter of fact those cases for which 
» = 1 and 2 are omitted in the table to save space. A missing » has » = 1 or 2 according 


| as (p? — 1)/8 is odd or even. This table is an extension of the large table of » for p < 300000 


published by M. Kraitcuik, in his Recherches sur la Théorie des Nombres, v. 1, Paris, 1924, 


| p. 131-191. Tables of v, like so many other difficult number-theoretic functions, must be 


constructed one entry at a time, and each entry must stand on its own merits, without 
hope of a smoothness check. Thus a small table like the present one represents a considerable 
amount of careful computation. 
The writer has in his possession a photostat copy of a manuscript containing an exten- 
sion by Beeger of the table here reviewed; it covers the primes between 310000 and 320000. 
Arithmeticians will be glad to learn that Mr. Beeger not only survived the war but 
was able to continue his researches in the theory of numbers. 


D. H. L. 


280[F].—A. GLopEN, “Table des solutions de la congruence X4+1=0 
(mod ~) pour 2-105<p<3-10°,” Mathematica (Rumania), v. 21, 1945, 
p. 45-65. 


This table is an extension of two previous tables: 
A. J. C. CunnincuaM, Binomial Factorisations, v. 1, London, 1923, p. 1-95, and v. 2, 1924, 
p. 1-104, for » < 105; and S. Hoprenot, “Table des solutions de la congruence 
x‘ = — 1(mod N) pour 100000 < N < 200000,” Brussels, Librairie du Sphinx, 1935, 18 P.; 
see MTAC, RMT 48, v. 1, p. 6. 

In order that the congruence 
(1) 1 =0 (mod p) 


have a solution it is necessary and sufficient that the prime p be of the form 8n + 1. Ac- 
cordingly, the table is concerned with such primes only. These*primes are capable of 
being represented in essentially one way by the quadratic forms 
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If one solves the linear congruences 
x =uy(mod p) and z = ot (mod p) 
for u and 2, then the four solutions of the quartic congruence (1) are given by 
X = +u(v + 1). 


The table gives for each » = 8n + 1 between 200000 and 300000 two values of X 
which are less than p/2; the other two values of X are of course negatives of these modulo p. 
There is one curious exception at = 200569 where, for some reason, the second solution 
is given as 100317 instead of 100252. Since the product of all the roots of (1) is clearly 
congruent to unity modulo #, the product of any two roots is congruent to + 1 modulo p. 
This fact was used as a check of the table. 

At the end of the table is given a list of 25 numbers of the form x* + 1 and two of 
the form x* + 1 with their complete decompositions into primes made possible by this 
table. For example j 

16204 + 1 = 289889-23759009. 


The first factor follows from the entry 
b = 289889, X, = 1620, X2 = 15926 


of the main table. The second factor is clearly a prime having no factors < 300000. A more 
extensive factor table of 54 entries (including only 8 of the above 27) based presumably 
upon a corresponding table of solutions of (1) up to p < 500000, appeared in Mathesis, 
v. 55, 1945, p. 81-82. 

It is perhaps worth pointing out that the main table is useful as a list of primes of the 
form 8 + 1 even if one is not concerned with the quartic congruence (1). 


D. H. L. 


281(F].—V. Tuésautt, Les Récréations Mathématiques (Parmi les nombres 
curieux), Supplément 4 Mathesis, Gembloux, p. 1-80, v. 54, no. 10, 
Aug. 1943; p. 81-119, no. 11, Feb. 1945. 16.325.2 cm. 


Here are numerous independent tables, not all published for the first time, including 
the following: 


1. With the digits 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, taken once, form two numbers of three digits 
whose sum is equal to numbers of (a) three digits chosen once; (b) four digits chosen once. 
On p. 8, for (a) a table of 21 solutions is given, and for (b) 48. For example, 286 + 173 = 459; 
789 + 246 = 1035. For (a) see L’@ducation Math., v. 23, 1920, p. 34-35. 

2. With the digits 1, 2, 3, 4, 5, 6, 7, 8, 9, taken once, form an integer which is a square. 
On p. 11, 30 solutions are given from 11826? = 139854276, to 30384? = 923187456; 29 of 
these were published earlier in Mathesis, v. 51, 1937, p. 279-281. On p. 12 are 87 squares 
formed from the digits 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, taken once, from 32043? = 1026753849, 
to 99066? = 9814072356. 

3. P. 62-68 gives ali 1-digit, 2-digit, 3-digit, and 4-digit square endings for the base 10, 
together with an inverse table showing the forms of those numbers whose squares end in 
a specified 1-digit, 2-digit or 3-digit number. The author has published this before in 
Mathesis, v. 48 Suppl. Oct. 1934, 22 p. A slightly more complete table is due to ScHapy, 
J. f.d. reine u. angew. Math., v. 84, 1877, p. 86-88. 


4. “‘Table des carrés des nombres entiers de 1 4 1000 dans les systémes de numération 
de bases B = 2, 3, 4, 5, 6, 7, 8, 9, 11, 12,” p. 89-119. 


These four tables are of interest to “digitologers,”’ a sect of which the author is a leading 
exponent. Table 3, however, is of practical use in the rapid identification of non-squares. 


D. H. L. 
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282[I]). , “Formulas for direct and inverse interpola- 
tion of a complex function tabulated along equidistant circular arcs,” 
J. Math. Phys., v. 24, Nov. 1945, p. 141-143. 17.4X26.4 cm. 


The purpose of this paper is to facilitate direct as well as inverse interpolation of an 
analytic function f(z) whose numerical values are known in a sequence of points z, which 
are equally spaced along a circle in the complex plane. In terms of the new variable 
P = (2 — x)/(z: — 2) the problem reduces to the case when the functional values are 
given at the points P, defined by 


= +0”) 1), @ = e*). 


v=p—n 
be the n-point Lagrange-Hermite interpolation formula based on our points P,. Setting 
Pr = >, + igqr (r = 0, 1, 2,---), the paper furnishes explicit expressions, in terms of 
pr, Gr and the angle 0, of the real and imaginary parts of the coefficients L,™(P), for 
n = 3, 4, and 5. It is also shown how these expressions may be used to apply also in the 
present situation the author’s recent formula for inverse interpolation (see “‘A new formula 
for inverse interpoiation,” Amer. Math. So., Bull., v. 50, 1944, p. 513-516). 


I. J. SCHOENBERG 
Univ. of Pennsylvania 


283(I].—HERBERT E. Sauzer, “Inverse interpolation for eight-, nine-, ten-, 
and eleven-point direct interpolation,” J. Math. Phys., v. 24, May 1945, 
p. 106-108. 17.4 26.4 cm. 


This paper extends a formula which was first published by the author as “A new 
formula for inverse interpolation,” Amer. Math. So., Bull., v. 50, 1944, p. 513-516. The 
expression in question is the inversion of the Lagrange interpolation formula. 

If we write f. = f(ao + xh), where h is the tabular interval, and if p is desired when 
we are given f,, then p can be expressed in terms of the tabulated values, fo, fi, fo, --- 
fu, f-2, etc. by means of the formula: 


m=1 da" 


in which L;™(qa) is a polynomial of (mn — 1) degree in a. 
The author reduced this expression to the form 


p=r— rs + — t) + + Sst — u)+-- 


where the quantities 1, s, ¢, u, etc. are expressed in terms of the tabulated values of the 
function. 

In the present paper the explicit values of these variables are given for the case where 
eight-, nine-, ten-, or eleven-point direct interpolation is used. 

It is difficult to imagine when an equation of this extent may be required in practical 
work, but the author remarks: “‘its full use can provide unusual accuracy in solving equa- 
tions (both real and complex) up to the tenth degree when the values of the polynomials 
are tabulated near the root at equal intervals.” 

H. T. D. 


284[I]—HeERBErT E. Sauzer, ‘‘Table of coefficients for double quadrature 
without differences, for integrating second order differential equations,” 
J. Math. Phys., v. 24, Nov. 1945, p. 135-140. 17.4X26.4 cm. 
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A function f(x) may be expressed in terms of its second derivative as 
(1) flee) = flea) + (x — + der 
zo 
If numerical values of f’’(x) are given in the m points 


Xp=x+ ph where = [n/2]), 


we may approximate f’’(x) by the m-point Lagrange formula. A double integration of the 
interpolating polynomial leads from (1) to the approximation formula 


@) + BR) = fle) + Con) + + + Ry 
i=u—n 


where C;,,(p) are polynomials in p of degree n + 1, given by 


the prime on the product sign indicating that the factor p2 — i is omitted. The paper 
furnishes for m = 3(1)11, the numerical (rational) values of the coefficients (3) for 
b = (u — nm + 1)(1)u, reduced to their least common denominator. 
The author discusses briefly the application of the formula (2) to the numerical inte- 
gration of a differential equation of the type 
+ o(x, y) = 0, 
to which any second order linear differential equation may be reduced by a proper change 
of unknown function. The solution once “‘started,”’ the familiar procedure of the progressive 
refinement of a first approximation to a new value of y applies, by working back and forth 
between the columns of values of y and y’’. At each step, however, the value of hf’(xo) 
must also be determined by (2) for any single value of p (#0) for which f(xo + ph) is 
already known. Thus a column of values of hy’ is also determined. The author concludes 
by pointing out further useful applications to the integration of partial differential equa- 
tions of the type 
= Ur + x, t). 

I. J. SCHOENBERG 


285(K].—@yvinp Burrau, “Middelfejlen som Usikkerhedsmaal,” Mate- 
matisk Tidsskrift, 1943B, p. 9-16. 15.124 cm. 


Suppose «1, X2, +++, X» are values of x in a sample of m items from a normal population 
with mean a and standard deviation ¢. W. S. Gossett, under the pseudonym ‘Student,’ 
in 1908 conjectured that the sampling distribution of 


2—a 
y (xi — 2)? 
was given by? 
(1) f(z)dz = C(1 + 2*) dz 


where C is a constant, such that | f(z)dz = 1. This conjecture was verified in 1925 by 
R. A. FIsHER.’ Gossett originally tabulated values of f. f(v)dv for which z = .1(.1)3, and 
n = 4(1)10. Fisher* has tabulated values of t = eV (n — 1) for which 


+2y(n—1 
= .1(.1).9, .95, .98, .99, 
(-1).9, .95, .98, .99 


and for = 1(1)30. 
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Mr. Burrau tabulates values of the integral 
(n—-1) 

2 d 

(2) 

for which z\¥(m — 1) = 0(.2)5, and for m = 4(2)12 and «. The case for m = = is, of course, 

1 
equivalent to finding the values of the Gaussian integral = f é edz for t = 0(.2)5. 


He then observes that since the mean value of 2’, i.e. E(z*), is 1/(m — 3), thus making 
E{(n — 1)2?] = 1, one might consider the tabulation of the integral 


3 v)dv 
(3) 
for ov (n — 3) = 0(.2)5, and for m = 4(2)12 and ©, rather than the analogous tabulation 
of (2), with the hope of obtaining values which change relatively little as m changes. The 
tabulation of (3) is made and the resulting values are considerably more constant with 
respect to changes in m than the analogous tabulations of (2). For example, the values of 


(2) and (3) for — 1) = 2 and (nm — 3) = 2are: 


n Value of (2) for sV¥(m —1) =2 Value of (3) for s¥(m —3) =2 
4 0.861 0.959 

6 0.898 0.950 

8 0.914 0.950 

10 0.923 0.950 

12 0.929 0.954 

oo 0.954 0.954 


This indicates that the “significance level’’ corresponding to + v(m — 3) is approxi- 
mately 95 per cent for ail values of nm = 4. The author does not discuss the problem for 
n <4. 

S. S. W. 


1 W. S. Gossett (Student), (a) ““The probable error of a mean,” Biometrika, v. 6, 1908, 
p. 1-25; (b) Collected Papers, Cambridge, 1942, p. 11-34. 

2 The z being used here is that originally used by “Student,” and is being used through- 
out this review. Mr. Burrau, although actually writing about “Student’s” problem, unfor- 
tunately denoted z/n* (in “‘Student’s”’ sense) by the letter z. 

oan” Fisuer, “Applications of ‘Student’s’ distribution,” Metron, v. 5, no. 3, 1925, 
p. 1 

4R. A. FisHer, ibid.; the table has also been published in Fisher’s Statistical Methods 

for Research Workers, eighth ed., Edinburgh, Oliver & Boyd, 1941, p. 167. 


286[K, M].—A. H. R. Grimsey, “On the accumulation of chance effects 
and the Gaussian frequency distribution,’ Phil. Mag., s. 7, v. 36, Apr. 
1945, p. 295. 16.8 25.3 cm. 


There are Tables of Pn, Qn, Un, for 
being both given for P,, where 


2 sinv\" 


Q, are the values to which P, approximates. 

Mr. Grimsey’s communication was inspired by misprints in L. SILBERSTEIN, “The 
accumulation of chance effects and the Gaussian frequency distribution,” Phil. Mag., 
s. 7, v. 35, 1944, p. 395-404, and in particular by his values of P, for m = 1(1)6, accom- 
panied by the statement “‘so far as we know, these integrals have not been previously 
determined,” and by the suggestion that mechanical quadrature would give the values for 
larger n. Mr. Grimsey refers to detailed discussion bearing on this question in J. EDWARDs, 
A Treatise on the Integrai Calculus, v. 2, London, 1922, p. 204-212, where the values of P,, 


ll 


[1(1)12; 8D], exact rational and decimal values 
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are incidentally given for = 1(1)7, that for Ps = 11/20 or 88/160, occurring in place of 
83/160 given by Silberstein on p. 399. 


R.C. A. 


287[L].—Mur.an S. Corrincton, “Tables of Bessel function J,(1000),”’ 
J. Math. Phys., v. 24, Nov. 1945, p. 144-147. 17.3 26.3 cm. Compare 
MTE 74, v. 2, p. 47-48; the Meissel errors are again noted here. 


Quotations: A table of Bessel functions for J,(1000) is given for n = [935(1)1035; 8D], 
also a graph of J,(1000) for 930 < m < 1040. All computations were made to 11D and 
later rounded off to 8D. 

When a frequency-modulated radio wave is analyzed into a spectrum, the carrier and 
side-frequency amplitudes are proportional to Bessel functions of the first kind, with order 
equal to the number of the side frequency and argument equal to the modulation index, 
in accord with the following equation: 


sin (wet + (D/u) sin 2xut) = Jn(D/u) sin (we + 2nnu)t 


where 

Wo = carrier frequency, radians per second, 
t = time in seconds, 

D = deviation in cycles per second, 

» = audio frequency in cycles per second. 


According to the present standards for commercial broadcasting, the deviation can be as 
much as 75,000 cycles per second and the audio frequency can be as low as 30 cycles per 
second, so the maximum value of the modulation index, D/p, is 2500. The table in this 
article was computed to determine the relative side frequency amplitudes near the edge of 
the band for a modulation index of 1000. The graph shows that when the side frequency 
number n exceeds 1000, the amplitudes decrease very rapidly towards zero. 


288[L].—J. Cossar & A. Erpéty1, Dictionary of Laplace Transforms. 
Admiralty Computing Service, Department of Scientific Research and 
Experiment, London; Part 3A, no. SRE/ACS. 102, Nov. 1945, 29 leaves. 
20.2 X33 cm. Mimeographed on one side of each leaf. These publica- 
tions are available only to certain Government agencies and activities. 
This penultimate part of the Dictionary is section VII. Section VI (Parts 2A, 1944; 
2B, 1945) was a classification according to f(t), the Laplace transform of f(t) being defined 
by the relation 

see M7AC, v. 1, p. 424-425. Part 3A is a Table of inverse Laplace Transforms, that is, a 
classification according to ¢(p) for the Elementary Functions, Gamma and Related Func- 


tions; and Functions defined by Integrals—as in Part 2A. Presumably Part 3B will classify 
according to ¢(p) the other functions listed on p. 425. 


R. C. A. 


289(L].—Rurvus Isaacs, ‘Airfoil theory for flow of variable velocity,” 
J. Aeronautical Sci., v. 12, 1945, p. 117. 20.5X28.6 cm. 


It is stated that Theodorsen’s tabulation! of the function 
C(x) = + iHo®(x)] 


is incorrect for x < .1 and so a new table is given for the real and imaginary parts of 
C(x) = F + 1G, for x = 0(.02).1(.05)4. There are 3D for F, 4D for G. 
H. B. 
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1 THEODORE THEODORSEN, “General theory of aerodynamic instability and the mecha- 
nism of flutter,” NACA, Techn. Reports, no. 496, 1934, p. 418, 426; see MTA C, v. 1, p. 404. 


In MemoriaM Note: This review was included in a letter written to me on Janua 
1946, the day on which Harry Bateman left Pasadena, California, by train c= York? 
where he was to have received notable honors. On the following morning his rare spirit 
took flight, as a result of a coronary thrombosis. The loss to scholarship, to our Committee. 
and to MTAC, of this extraordinarily able, profoundly learned, and ever helpful intellec- 
. B. was born in Manchester, England, in 1882. He graduated A.B., as senior 
at Trinity College, Cambridge, in 1903, and was a Smith’s Prizeman in 1905. During 19082 
1911 he was a fellow of Trinity and received his A.M. in 1906. He was a student at Géttingen 
and Paris 1905-06, a lecturer at the University of Liverpool 1906-07, and a reader of 
mathematics and physics at the University of Manchester 1907-10. Coming to the United 
States in 1910 with 50 publications already to his credit (they included the notable report 
on Integral Equations), he was for two ong a lecturer at Bryn Mawr College. Duri 
1912-15 he was a Johnston scholar at The Johns Hopkins University (Ph.D. 1913), an 
was lecturer there from 1915 until 1917, when he was appointed professor of mathematics, 
physics, and aeronautics at the California Inst. Tech. He was elected a fellow of the American 
Philosophical Society in 1924, of the Royal Society of London in 1928, and of the National 
Academy of Sciences in 1930, and a vice-president (1935) and Gibbs lecturer (1943) of the 
American Mathematical Society. His published books were The Mathematical Analysis of 
Electrical and Optical Wave- Motion on the Basis of Maxwell's Equations, Cambridge, 1915; 
Differential Equations, London, 1918; Partial Differential Equations of Mathematical Physics, 
Cambridge, 1932; new edition, New York, 1944. He edited the lectures of H. A. LorENTz, 
Problems of Modern Physics, Boston, 1927. He was also a collaborator on three published 
reports of the National Academy of Sciences: (a) Committee on Electrodynamics of Movin 
Media, 1922; (b) Committee on Hydrodynamics, 1932; (c) Committee on Numerica 
Integration. His report on tables of Bessel functions was published in MTAC, no. 7, 1944 
and our Committee had hoped soon to receive his manuscripts of reports on tables of Higher 
Functions and Integrals, and of a great work on Definite Integrals, replacing the antiquated 
volume of BIERENS DE HAAN; compare MTAC, v. 1, p. 322. About 1900 H. B. participated 
in a chess tournament between England and America; and prior to 1911 he beat Steinmetz 
in a tournament. For other details see Who’s Who 1945, London; Who's Who in America 
1944-45 ; Amer. Men of Science, seventh ed., 1944; Who's Who in Engineering, 1941; Who's 
Who in Aviation. . .1942-43; Pan Pacific Who’s Who,. . . 1940-41 edition, Honolulu, Hawaii 
1941; Leaders in Education, second ed., 1941; C. M. Neate, The Senior Wranglers of the 
University of Cambridge, Bury St. Edmunds, 1907; “Poggendorff,” v. 5, 1925, and v. 6, 1936. 
R. C. A. 


290[L]—_NYMTP, Table of Characteristic Values of Mathieu's Differential 
Equation. Prepared for the Applied Mathematics Panel of the National 
Defense Research Committee. No. AMP Report 165.1 R. New York, 
September, 1945. xxiv, 39 p. 20X28 cm. Produced by photo-offset 
from typescript. These tables are available only to certain Government 
agencies and activities. 


The aim of the volume is to provide an 8D table of the characteristic values of Mathieu’s 
differential equation, for orders up to 15, interpolable in the parameter s (defined below), 
within the range given. 

An account of previous tables of these values has been given by the reviewer in MTAC, 
v. 1, p. 409-419. 

Mathieu’s differential equation is written in the form 


(1) d*y/d? + (b — scos*t)y = 0. 


For a prescribed value of s this equation has a solution admitting the period 2x in ¢ only 
when 6 has one of a discrete but infinite sequence of “characteristic values.” The periodic 
solutions fall naturally into two classes, according as they are even or odd in ¢, and the 
corresponding characteristic values are denoted respectively by be, and bo,, the suffix 
denoting the order in the sequence arranged according to increasing numerical magnitude. 
The ranges included are r = 0(1)15, 0 < s < 100; the interval varies, being chosen so 
that the maximum attainable accuracy is obtained by the use of modified second differ- 
ences, and these differences are provided. In detail 
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be, bo, 
0 0(.2)20(.5)37(1)100 1(1)4 0(.5)100 
1 0(.2)15(.5)50(1) 100 5(1)7 0(1)100 
2 0(.2)40(.5)77(1)100 8,9 0(2)100 
3 0(.2)13(.5)56(1) 100 10(1)15 0(4)100 and 10(20)90 
4,5,6 0(.5)100 
7,8 0(1)100 
9, 10 0(2)100 
11(1)15 0(4)100 and 10(20)90 


A table of the Everett interpolation coefficients, E2(p) and F2(p), for p = [0(.001)1; 6D], 
is given. 

A graph showing the characteristic values as functions of s is reproduced on p. xxiv of 
the introduction. 

It should suffice to say that the NYMTP has carried out the task assigned to it in the 
efficient manner we have come to expect from this Project. It is clear from the introduction 
that considerable advances have been made in the technique for the routine computation 
of these characteristic values, although the new methods do not become really powerful 
until several decimals of. the characteristic value have been determined. In this respect 
the computer derives, from the small interval, advantages in extrapolation, commensurate 
with those in interpolation which he confers on the user. Despite some adverse criticism 
to be given below, the Table is to be welcomed, since any table is infinitely better than no 
table, and because this seems to be the first table of any magnitude which is easily inter- 
polable in the parameter, and in which provision for interpolation is made. It is also funda- 
mental, since a knowledge of the characteristic values is necessary for the calculation of 
the Fourier coefficients, and hence of the periodic functions themselves. 

On the evidence available to us, however, it does appear that the NYMTP was upon 
some points ill advised (or instructed). Through ignorance or perversity the judgment of 
Matuieu, HEINE, WHITTAKER, GOLDSTEIN, INCE, and nearly every other worker in the 
field, has been ignored or disregarded by choosing the form (1) of the differential equations, 
a form which is inconvenient in so many ways both theoretically and practically, and which 
seems to have only one trivial circumstance in its favor. Indeed, the explanation of the 
methods of computation and the formulae given in the introduction, and the computations 
themselves, are all relative to the equation 


(2) dy /dt? + (a — 26 cos 2t)y = 0 
in which, of course, the parameters are related to those in (1) by 
(3) a=b—4ts, 


The superiority of (2) over (1) as a canonical form seems to be almost universally rec- 
ognized. 

Fundamentally this rests upon the simplicity and symmetry of the relations between 
solutions for positive and negative values of @ (or s). For this, reference may be made to 
MTA C, v. 2, p. 1-11. Changing the sign of @ in (2) leaves a unchanged, and is equivalent 
te replacing ¢ by + 34 + t. Had @ and not b been tabulated, the table would have served 
equally well for negative and positive values of s, but as it is, it is necessary, in order to 
derive the characteristic number corresponding to a negative value of s, to subtract |s| 
from the value corresponding to |s|. 

This point can be seen more immediately and strikingly if one compares the elegant 
symmetry of the graph of characteristic values against @ (or some equivalent parameter) 
as given by IncE,! Srrutt,? Lusxin & SToker,* McLAcuian,‘ and others, with the dia- 
gram on p. xxiv of the present volume, imagined carried into the region of negative s. 

Moreover, we are so far given only characteristic values. If the Fourier coefficients of 
the periodic solutions are needed—essential before the functions themselves can be calcu- 
lated—then all the formulae are most conveniently expressed (as in the introduction !) in 
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terms of a. All the indications go to show that NYMTP did, in fact, compute a, and then 
converted to b—and that in computing the tables of the corresponding Fourier coefficients 
(which are promised) it will be a and not 6 that NYMTP will actually use. 

It is fervently to be desired that information equivalent to that given in these tables 

may before long be made more generally available, but we urge, with all the force at our 
command, that before this is done the tables should be reconverted so as to give a and not b. 
In one respect, perhaps, equation (1) gives a lead that was no longer followed. It indi- 

cates more clearly than does (2) that the parameter s is generally connected with the square 
of some significant physical quantity. This is an argument for the retention of the param- 
eter c (s = c*) of some earlier tables.5 

I, W. G. BICKLEY 

Imperial College of Science and Technology, London 

of 1E. L. Ince, “Researches into the characteristic numbers of the Mathieu equation,” 
R. So. Edinburgh, Proc., v. 46, 1925, p. 28. 

he 2M. J. O. Srrutt, "Lamésche, Mathieusche, und verwandte Funktionen in Physik und 

on Technik, Berlin, Springer, 1932, p. 24. 

3S. Luskin & } “Stability of columns and strings under periodically 
oe varying forces,” Quarterly Appl. Math., v. 1, 1943, p. 221. 
ful 4N. W. McLacaian, “Computation of the solution of Mathieu’s equation,” Phil. Mag., 
act s. 7, v. 36, 1945, p. 404. 
te bd STRATTON, Morse, Cau & Hutner, Elliptic Cylinder and Spheroidal Wave Functions, 
om New York, Wiley, 1941; compare MTAC, v. 1, p. 157f, 414. 
no 
a. 291[L]—NYMTP, Tables of Associated Legendre Functions. New York, 
la- Columbia University Press, 1945, xlvi, 303 p. 18.626.5 cm. Repro- 

of duced by a photo-offset process. $5.00. 

Another important volume has been added to the impressive list of publications of 
on the NYMTP. This work, devoted to the associated Legendre functions of first and second 
a kind, is a very useful addition to the already extensive number of tables which have ap- 

peared in this area of functional computation as described earlier in MTAC, v. 1, p. 116- 
~ 119. The present table ‘‘was begun in 1940, to meet urgent needs for a table to about six 
rn significant figures at intervals of 0.1. At this interval in the argument, interpolation in 

these functions is not everywhere satisfactory; a very considerable amount of further work 
- would be required, in order to make all these functions reasonably interpolable over the 
entire range covered.’’ The Project has on its agenda the subtabulation of the functions to 

a more satisfactory interval of the argument. 

By the associated Legendre functions of mth order and mth degree are meant the 

functions P,."(x) and Q,(x), the former called the first kind and the latter the second kind. 

In the range x? = 1, it is customary to define these functions as follows: 

P,™(x) = (1 » amx) =(1- 
een 
to where P,,(x) is the nth Legendrean polynomial, and Q,(x) the mth Legendrean function of 
ent second kind, defined as follows: 
ved 
1 
3(n — 1) 5(m — 2) 
lia- It should be observed, however, that H. TALLevist, in his tables, defines Q,(x) as twice 
the value of the one given. 
; of In the range x? > 1, it is customary to change the multipliers of the derivatives of 
cu- P,(x) and Q,(x) to (x? — 1)!" and to replace in the definition of Q,(x), the function 
} in log [(1 + x)/(1 — x)] by log [(x + 1)/(« — 1)}. 


© 
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In all previous computations of the Legendrean functions, the range of the independent 
variable has been restricted to the interval — 1 =x = + 1. The present work thus makes 
a fundamental contribution, since most of the tables are computed over ranges which lie 
outside the unit interval. 

The present volume contains 19 tables which we may describe as follows: 

In the first ten tables the values of m and m are taken over the range m = 1(1)10; 
m = 0(1)4, m =n, and the computation is carried to 6S. I. P,™(cos 6), @ = 0(1°)90°. 
Il. dP,™(cos 6)/de, @ = 0(1°)90°. I. P,™(x), x = 1(.1)10. IV. dP,™(x)/dx, x = 1(.1)10. 
V. (— 1)™Q,™(x), = 1(.1)10. VIL (— 1)™*dQ,™(x)/dx, x = 1(.1)10. VIL. i-*P,™(ix), 
x = 0(.1)10. VII. i-*dP,™(ix)/dx, x = 0(.1)10. IX. i**™4Q,™(ix), x = 0(.1)10. X. 
(ie) x = 0(.1)10. 

Of the remaining nine tables in the volume, the first four give the values of the func- 
tions and their first derivatives at fractional half orders, that is to say, for the subscript 
m+ 4, = — 1(1) + 4. In these tables m = 0(1)4. The tables may then be described as 
follows: 

XI. Priy(x), x = (1(.1)10; 4S-6S]. dPt,4(x)/dx, x = [1(.1)10; mainly to 6S]. 
XMM. (— a(x), = [1(.1)10; 6S]. XIV. (— x = [1(.1)10; 6S]. 

The remaining five tables are supplementary and are contained in 10 pages. They 
may be described as follows: 

XV. Exact values of P,™(x) and dP,™(x)/dx, m = 0, 2, 4; m = 4 or 5(1)10, x = 1(.1)2. 


XVI. Exact values of P,™(x)/(x? — 1) and of [Srrc@|/e — 1), m =1,3,n=2 


or 3(1)10, x = 1(.1)2. XVII. Values of (x? — 1)#P%,4(x) and of (x? — 1)*dPiy4(x)/dx, 
m= 1, 3, m = — 1(1)4, x = [1(.1)2; 4D-7D]. XVIII. Values of Legendre’s normalizing 
factor, N,™ = {2(n + m)!/[(2n + 1)(n — m)!]}4, n = 0(1)10, m=n, 10S. XIX. First 
eleven coefficients in the expansion of CF(m + 4; } — m; n + 3; — 4), where C is equal 
to + m + + HI}. 

Except in the case of P,(cos @) and its derivative, which were computed directly from 
their expansions in terms of sin ~@ and cos ~@, the other tables were evaluated by means of 
standard recurrence formulae. In the final checking of the computation the following 
formula was used: 

— (— 4)"T [k(n + m + + m + 2)] 

(x? — 1)T — m + — m + 2)] 
In the case of P,,™(cos 0), when the values were checked against those published by Tall- 
qvist, no errors were discovered. 

In addition to the tables, a number of formulae are developed which were found useful 
in the computation of preliminary values. A bibliography of 47 references is also given, 


which provides an adequate survey of the present status of computation in this active field. 


d 


292[L, M].—P1ERo GiorGio Borponl, conical sound source,” Acousti- 
cal Soc. Amer., J., v. 17, Oct. 1945, p. 125. 


There is a table of G,(z), for m = 0(1)6, z = [0(1)5; 3D-4D]. 


293[L, M].—Great Britain, NAUTICAL ALMANAC OFFICE, Tables of the Inte- 
grals C(t) cos (2tk'u —ku*)du, S(t) =k! fore sin (2tk'u —ku*)du, 
Department of Scientific Research and Experiment, Admiralty Com- 
puting Service, September 1945. No. SRE/ACS. 90. 11 p. 20.1 x 30.8 
cm. This publication is available only to certain Government agencies 
and activities. 
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Tables to 3D are given for C, S, and (C? + 5S*)4, for the following values of & and ¢: 


k=i3, #=— 10, — 3(.1) + 3, 10; 
k=}, t= —8(1) — 5(.5) — 3(.1) + 8; 
k=1, t= —6, — 55) — 3(.1) + 6; 
k =2,4,t = — 5(.5) — 3(.1) + 7; 
k=8 t= —5(.5) — 3(.1) +9. 


“Graphical accuracy only was required and the method of computation was such that the 
last of the three decimals tabulated may be in error by a few units.”” There is a graph of 
(S+ C)tfor —S<t< +6. 

The differential equations satisfied by S(#) and C(é) are 


dS/dt = 2tC — S/ki, dC/dt = 1 — 2tS — C/ki, 


for the rangé of values — 3 < ¢ < + 3; S(t) and C(t) were obtained by numerical inte- 
gration of these equations. For larger values of t, S(f) and C(t) were computed from certain 
appropriate asymptotic trigonometric series. 

R. C. A. 


294[M].—W. HEISENBERG, ‘“‘Die beobachtbaren Gréssen in der Theorie der 
Elementarteilchen. II.,’’ Z. f. Physik, v. 120, 1943, p. 702. 


There is a short table of 


pi = 1, po = 3, ps = 13, Pe = 79, D = GA1, Ps = 6579, etc. It is, of course, remarkable 
that these integrals are integers. 
H. B. 


295[U].—HAmBurRG, DEUTSCHE SEEWARTE, publication no. 2154, F-Tafel. 
Tafel zur vereinfachten Berechnung von Héhenstandlinien. 3 Auflage. 
Hamburg, August, 1941. xxiii, 88 p. 19.6 X 29.2 cm. In the third 
edition there were extensions and corrections of the introductory mate- 
rial, and of 8 of the 11 tables. 


The method and principal table of this volume are similar in many respects to those of 
H. O. 208 (Dretsonstox, see MTAC, v. 1, p. 79f). The astronomical triangle is divided 
into two right spherical triangles by a perpendicular from the zenith upon the hour circle 
of the star; U is the co-declination of the foot of the perpendicular, and V is log cos B, 
where B is the angle subtended at the zenith by U. By Napier’s rules, 


tan U = costcot L 
and 
sin B = sin ¢ cos L, 


where ¢, L, and d are the local hour angle, latitude and declination respectively. By apply- 
ing another of Napier’s rules to the right triangle of which the star is one vertex, the alti- 
tude, h, may be found by 
sin h = cos B sin (d + U) 
or 
log sin kh = V + log sin (d + U). 


For the determination of azimuth, Z, two more auxiliary quantities are introduced, 
P which is the great circle distance from the star to the east- or west-point of the horizon, 
and Gr. 6 which is the declination of the intersection of the hour circle of the star with the 
prime vertical. Thus, sin¢cosd = cos P and sin Z = cos P sech. Also, tan Gr. 5 = tan L cost. 
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In Table F I, with vertical argument, latitude 0(1°)70°, and horizontal argument, 
local hour angle 0(4™)6", three values per page, there are tabulated four quantities, U to 
the nearest 0’.1, V to 5D, Gr. 5 and P, each to the nearest 0°1. In the second part of 
Table F I, the vertical argument is latitude, 70°(1°)90°, and the horizontal argument is 
local hour angle 0(4™)6", nine values per page, three in each horizontal section. 

At the bottom of the vertical columns in Table F I are azimuths; entering the left 
hand column with altitude as argument, and moving across the pages horizontally until one 
finds under P, the value already copied out, one can drop to the bottom of the column 
and read off the azimuth angle. Since the tabulated values of the azimuth angle go up to 
90° only, it is necessary to have another device to determine the quadrant. When the 
hour angle is greater than 6", the azimuth is measured from the elevated pole; when the 
local hour angle is less than 65 and L and d are of opposite name, the azimuth is measured 
from the depressed pole. If the local hour angle is less than 6" and L and d are of the same 
name, the azimuth is measured from the elevated or depressed pole according as the decli- 
nation is greater or less than the quantity Gr. 6. : 

In case the altitude is great, or the azimuth near 90°, the value of the azimuth may be 
poorly determined by the use of Table F I. In such a case, it will be noted that the value 
of P lies below a dotted line running across the page. One must then use instead Table 
F XI, which gives P to the nearest minute of arc and the variation in P corresponding to 
1’ change in d or h. 

Table F II is a table of log sin x, x = [0(0’-1)6°(1’)90°; 5D], with generous tables of 
proportional parts. 

Tables F III and F IV represent the principal advantages this volume possesses over 
other similar tables; they permit one to determine the corrections (to the nearest 0/1) to 
the computed altitude corresponding to slight changes in time (up to 2™ by 10° steps) or lati- 
tude (up to 30’ by 1’ steps) respectively. In both cases, one can interpolate very easily by a 
shift of the decimal point. Table F III is a well-designed triple-entry table occupying only 
five pages; one starts down the column at the left headed by the value nearest the assumed 
latitude, stops at the value nearest the computed azimuth and moves to the right to the 
column headed by the number of seconds change in time. Table F IV is a small double- 
entry table on a single page; the vertical argument is azimuth 0(5°)20°(2°)90°, and the 
horizontal argument is change in latitude, 1’(1’)10’(10’)30’. These two tables allow one to 
work either with an assumed position or with a dead reckoning position. 

Table F V is for changing time into angular measure and conversely. Table F VI gives 
the corrections for refraction, semi-diameter and parallax to be applied to the altitude 
(3°-90°) of the lower- or upper-limb of the moon; there is a supplementary table for height 
of eye. Table F VII gives the combined correction for refraction and height of eye (0-30 
meters) to be applied to the altitudes (3°-90°) of fixed stars or planets. Table F VIII yields 
the correction for refraction, semi-diameter and height of eye (0-30 meters) to be applied 
to altitudes (3°-90°) of the sun’s lower limb; there are also two auxiliary tables to provide 
corrections to the altitudes to take care of the varying semi-diameter of the sun through the 
year, and for the case where the sun’s upper limb was observed. The latter takes only a 
very small amount of space and would seem to be quite worthwhile. Tables F IX and F X 
provide similar corrections for use with the bubble sextant. 

The tables are well printed on a good grade of paper. In a number of cases, the rules 
needed to make decisions as to quadrants, etc. are printed on each page. As for the accu- 
racy of the tabulated values, only a few rounding off errors of a unit in the last place were 
discovered in a brief examination. 


CHARLES H. SMILEY 
Brown University 


296[U].—HyproGrapPuic Orrice, Publication no. 351, v. 7, Secret. Celestial 
‘Air Navigation Table. Volume 7, Latitudes 60°-70°, Tokyo, November 
1942. iii, 320 p. 18 & 25 cm. 
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This is the seventh volume of a series originally intended to include only six and to 
cover latitudes 60°S to 60°N. Since its general design and content are similar to those of 
the first six volumes which have already been reviewed (RMT 273), it will suffice to point 
out the only change which has been made. The auxiliary table for changing time measure 
into angular measure has been replaced by another which has been pasted over the old one. 
Only the degrees were given in the main body of the old table and one had to move to the 
right hand edge of the page to pick up the minutes associated with the degrees; in the new 
table, the degrees and the minutes are given together in each entry. 

It was earlier stated, MTAC, v. 2, p. 44, that v. 4 was the first of the series to be issued. 
Actually v. 3 appeared in August, 1940; v. 4 in May, 1941; v. 1 in September, 1941; v. 5 
in December, 1941; and v. 6 in February, 1942. A copy of v. 2 has not yet reached the 
reviewer's hands. 

CHARLES H. SMILEY 


OFFICE, Publication No. 9582. Secret. Altitude- 
Azimuth Tables for September 1945 (Sun, Moon, Planets and Fixed Stars) 
computed for Kisarazu, Naha, Iwo Jima and Kanoya and for Japanese 
Central Standard Time. Tokyo, Japanese Hydrographic Office, August 
1945, 42, vi, viii p., 18.2 X 25.7 cm. 


These tables were designed for celestial air navigation over a limited region in the 
western Pacific. All descriptive material, explanations, star names, etc., are in Japanese 
characters; tabular material is in Arabic numerals. One volume is prepared for each month 
and contains all the material needed for that month. For each 20 minutes throughout 
the 24 hours of Japanese Central Standard Time (time zone —9), altitudes and azimuths 
of the celestial bodies, most useful to navigators, as computed for four locations, namely, 
Kisarazu near Tokyo, Naha on Okinawa, Iwo Jima, and Kanoya near Kagoshima, are 
given to the nearest minute of arc and to the nearest degree respectively. For the sun, 
the altitude (> 4°) is given with refraction for 4000 meters elevation added. In the case 
of the moon, the altitude (also > 4°) with corrections added for refraction at 4000 meters 
and parallax, is given. For bright stars and planets, the altitudes (occasionally down to 2° 
and always with refraction for 4000 meters), are also given. 

Ten pages are used in presenting material for each of the four geographical positions 
listed above. Each page carries three vertical sections, each of which contains three double 
columns and includes data for a single day. The double column to the left gives the altitude 
and azimuth (say, 2956 050 which means 29° 56’ altitude, 50° azimuth) for a bright star 
until sunrise, then for the sun until sunset, and finally for another bright star until mid- 
night. The bright stars in this column are chosen so as to be useful in determining latitude. 

The second double column in the vertical section gives the altitude and azimuth of the 
moon when it is available for observation and for bright stars or planets for the remaining 
time. The third double column presents the altitude and azimuth of bright stars and planets. 

Since only the sun and moon are easily observed in the daytime, one or two double 
columns would ordinarily be blank, opposite the daylight hours as listed on the left side of 
the page. These spaces have been filled with data for bright stars and planets, well located 
for observation in the night hours; this material is printed in heavy black type and the time 
argument is printed on the right side of the page in heavy black type also. Thus there are 
available for observation, one or two celestial objects during the day and five or four during 
the night. The planets included in this volume are Venus and Saturn; the bright stars used 
are Capella, Vega, Alpheratz, Arcturus, Regulus, Altair, Procyon, Rigel, Spica, Sirius, 
Antares, Fomalhaut and Canopus. 

These tables differ considerably from common navigation tables in that much greater 
differences between the assumed and real positions of the observer are used here, and hence 
larger differences between computed and observed altitudes are encountered. In one of 
the examples worked out in the explanation in the back of the book, altitude intercepts of 
+2°07’ and —37’ are used in plotting two lines of position and determining a fix. The 
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inaccuracies introduced in this way would be too large for surface navigation, but may be 
tolerated in air navigation. 

In addition to the principal tables just described, there is one page on which the altitude 
of Polaris is given for each of the four geographical positions named and for each 20 minutes 
of the night for every fifth day of the month. There follows a one-page interpolation table, 
giving the corrections to be applied to tabulated altitudes. The vertical argument is tabular 
difference in altitude 0(4’)5°, and the horizontal argument is: time interval 0"5(0™5) 10™. 
There follow four pages of tables, one for each of the geographical positions giving correc- 
tions to the tabulated altitudes corresponding to other positions near the chosen ones. 
There are 10 such points near Tokyo, 8 near Iwo Jima, 9 near Kagoshima, and 6 near 
Okinawa. These positions vary as much as approximately 6° in latitude and 4° in longitude 
from the central position. It will be realized that these corrections will be most useful in 
the case of a body observed at low altitude. At the top of a column is the name of a geo- 
graphical location for which the altitude correction applies. The left-hand column of the 
page contains values of the argument Z (azimuth of the celestial body) 0(3°)180°; the 
right-hand column contains values of the argument Z 180°(3°)360°. No provision is made 
for determining the change in azimuth in going from the standard position to the nearby one. 

The remainder of the volume is devoted to explanations; an expert navigator may 
follow these, even though he does not read Japanese, since the numerals are Arabic. One 
notes that the page of tabular values reproduced for use in working the examples covers 
April 22-24, 1945, and that several improvements have been made in the arrangement of 
the tables since that time, particularly in the manner in which heavy black type is used to 
set apart night-time values printed opposite day-time hour arguments. 

Since much of the air travel over oceans will in the future be over a limited number of 
well-defined paths, it may be that the type of table described here may come into use over 
such paths. The advantage of dispensing with the need of an almanac and the rapidity with 
which a fix can be determined may be considered as outweighing in air navigation, the 
inaccuracies introduced by the necessarily greater differences between dead reckoning 
position and assumed position. There remains the question as to whether or not it will be 
economically profitable to prepare such tables, no matter how convenient and rapid they 
may be, if their useful life is only one month. 

CHARLES H. SMILEY 


298[U].—P. V. H. Weems, Line of Position Book, fourth ed. Annapolis, Md. 
Weems System of Navigation, 1943, viii, 56 p. 14.5 K 24.3 cm. $2.00. 


This book, in its first edition issued in 1927 by the U. S. Naval Institute, contained 
tables useful for latitudes 65° South to 65° North and was intended for surface navigation. 
The second edition appeared in 1928 and extended its field to include all latitudes; it also 
contained simplified procedures for aerial navigation. The third edition followed after an 
interval of 12 years. The author is now a Captain (retired) of the U. S. Navy and has long 
been well known for his contributions to the art and science of navigation. The tables in 
this volume are said to provide altitudes accurate to the nearest minute of arc without 
interpolation; greater accuracy can be had by interpolating. 

There are three principal tables, Table A, Table B, and Table A-Polar. Table A has 
as vertical argument, latitude 0(1°)65°, and as horizontal argument, local hour angle, 
0(1°)180° by integral degrees. The values, 0 to 90° of local hour angle are found at the 
tops of pages, 5° per page except the first page which has 6°, and the tabulated quantities 
are A and K; local hour angles 90° to 180° are found at the bottoms of pages, and A and 
(180° — K) are tabulated. In the notation of RMT 103, A is 105 log sec a and is given to 
one decimal place for local hour angles 0 to 20°, and 160° to 180° to the nearest integer for 
local hour angles 20° to 160°. K is the declination of the foot of the perpendicular from the 
zenith upon the hour circle through the celestial body and is given to the nearest tenth of a 
minute of arc throughout. 

Table B is a table of 10° log sec x for x = 0(1’)90° when the argument at the top and 
left of the page is used, or a table of 105 log csc x for the same interval when the argument 
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at the right and bottom of the page is used. Values less than 665 are given to one decimal, 
other values to the nearest integer. 

Tables A and B are credited to S. OGura, hydrographic engineer of the Japanese Navy. 

Table A-Polar is similar to Table A, except that latitudes 60° to 90° are covered and 
values of A less than 665 are given to one decimal, other values of A to the nearest integer. 

Weems suggests that the formula B(h) = A + B(K ~ d) be used for the determina- 
tion of the altitude and that the azimuth diagram of ArmistEAD Rust, USN Captain 
(retired), be used to determine the azimuth. This diagram is reproduced on a scale 1.5 
times the original on two pages, one of which is folded. A diagram enabling one to deter- 
mine the hour angle of a celestial body on the prime vertical takes care of the ambiguous 
case in the determination of the quadrant of the azimuth angle arising when the latitude 
is of the same name and greater than the declination. 

Auxiliary tables given in the volume are those providing corrections for observed alti- 
tudes obtained with ordinary sextants and with bubble sextants. A brief but adequate 
explanation of the use of the tables is given and several examples are worked out. 

In the preface to the third edition, the author says, “The Line of Position Book has 
stood the test of time. No errors have been noted in the Altitude Tables. Except for the 
method of H.O. 214, it remains the shortest, easiest method for working the line of posi- 
tion.”’ In this statement, the author fails to mention one of the strong points of his method, 
namely that it works in all latitudes and for all altitudes. All of the modern tables with 
which a sight can be reduced in appreciably less time than with the Line of Position Book 
are limited methods, usually limited in both latitude and altitude. 

Perhaps the outstanding criticism to be made of this book of tables is that local hour 
angle changes from page to page rather than latitude. There is a definite advantage in 
finding all of the data to be used at one time on a single page rather than scattered through 
the book. By the same argument, the polar table should be amalgamated with Table A. 
It may be noted here that these changes have been made by the author in his New Line of 
Position Tables, Annapolis, Md., 1944. 
CHARLES H. SMILEY 
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References have been made to Errata in RMT 286 (Silberstein); N 53 
(Furnas, Wrinch); 54 (McIntyre); Q 17 (?). 


79. Barlow's Tables of Squares, Cubes, Square Roots, Cube Roots, and Re- 
ciprocals, of all Integer Numbers up to 12500, edited by L. J. Comrir. 
Fourth ed., 1941. See MTAC, v. 1, p. 16-17. 


As a part of the checks applied in the transcription of these tables to punch-cards, the 
tabular values of f(m) = Vn, 107, Vn, and 1/n for » = 1000(1)12500 were inverted; 
i.e., f[f(n)] was found. Differential corrections accurate to a few additional decimal places 
were computed. Doubtful cases were further examined. All values were found to be correct, 
except two which need a correction of a unit in the last place: 

1) For n = 3195, V10n is tabulated as 178.745630. The correct value is 178.7456293, 
which rounds off to 178.745629. Also the corresponding differences must be changed to 
27975, 27971. 

2) For n = 4575, Vn is tabulated as 16.600851. The correct value is 16.6008516, 
which rounds off to 16.600852. Also the corrected differences, 1210, 1209. 


L. J. CUNNINGHAM 
Aberdeen Proving Ground 


EpiroriaL Note: These errors are also to be found in the third ed. 1930, and in the 1935 
pirated ed. 


80. James Burcess, “On the definite integral (2/4) ---," See 
MTAC, v. 1, p. 449. 


be 
ide 
tes 
lar : 
p=, 
ec- ; 
es. 
ear 
ide 
in 
20- 
the 
the : 
ide 
ne. 
ay 
ine 
ers 
of 
to 
of 
yer 
ith 
the 
ing 
be 
ey 
d. 
0. 
ed 
on. 
lso 
an 
ng 
in 
ut 
1as 
le, 
he 
ies 
nd 
to 
for 
the 
fa 
nd 
nt 


86 MATHEMATICAL TABLES—-ERRATA 


In trying to use Burgess’ Table of L as a check on certain computations we found 
several errors in this table. However, to decimals varying from 11 to 15, we checked 
Wrench’s corrections to Burgess’ table except for ¢ = 6. There we agreed with Burgess 
rather than Wrench. In order to settle this point once and for all we made three separate 
computations, one using the 15th convergent of a continued fraction for L/2t, the second 
using the 16th convergent, and the third using the asymptotic series for L. All three were 
sufficiently accurate to verify Burgess’ 15 decimal value and refute Wrench’s value. In 
addition we made use of a trick for increasing the accuracy obtainable from an asymptotic 
series to compute the value 

-98665 31092 31165 44447 for Latt = 6. 


J. BARKLEY RossEeR 
Epwarp M. Cook 
Joun A. BritTaIn 
Alleghany Ballistic Laboratory 
Cumberland, Md. 15. XI. 45 


I carefully recalculated L(6), this time using all 24 decimal places of F(x) and its 
reduced derivatives for x = 8.4 and 8.5 as tabulated by Sheppard. The final result is 
L(6) = .98665 31092 31165 44446 847 where the error in the 23rd place is less than one 
unit. This confirms the above result to 20D. In view of our agreement in all other cases 
this should conclude the checking of the L-table. 

J. W. WreENcH, Jr. 
4211 Second St., N.W. 
Washington 11, D. C. 20.X1.45 


81. OLINTHUS GILBERT GREGORY (1774-1841), Mathematics for Practical 
Men . . . second American ed. from the second London edition, Phila- 
delphia, Baird, 1863. In the supplementary tables is, “‘A Table of Circles, 
from which knowing the diameters, the areas, circumferences, and sides 
of equal squares are found.”’ At the conclusion of the tables Gregory 
states that they ‘“‘were computed with great care by the author’s esteemed 
friend, the late H. Goopwin, Esq. of Blackheath, a gentleman whose 
indefatigable perseverance and remarkable accuracy in reference to 
numerical computations cannot be too highly characterized. They are 
inserted here to supersede the necessity of consulting some erroneous 
tables of the areas, etc., of circles recently put into circulation.”’ 


These three tables are each for Diameter (D) = [1(.25)100.75; 8D], except that the 
Area for D = 1 is to 7D. With these same values of D I recalculated the whole of the 
table for circumferences (C) and found that of Goodwin’s 400 tabular values of C, 126 
are correct, one is too small by 3 units in the last figure, 22 are deficient by 2 final units, 
244 are too small by one unit in the last place, and one is too great by 3 final units. The 
following major errors were discovered in this table: 


D For Read 
37.75 118.59572267 118.59512267 
43.00 135.08348410 135.08848410 
48.00 150.79644797 150.79644737 
50.25 157.96503084 157.86503084 
81.75 256.82579942 256.82519943 
84.00 263.89378269 263.89378290 
96.25 302.37829269 302.37829291 


353 South Indiana St., 
Danville, Indiana 


STANLEY M. SHARTLE 
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Epitoriat Note: The posthumously published ““Goodwin” table, discussed by Mr. Shartle, 
appeared as Table II in the first London edition of Gregory’s work, London, 1825, p. 399- 
406; in the third London edition, 1848, edited, rev., and vat f by Henry Law, it is Table VII 
in the Appendix, p. 66-72. In each of these editions the seven sets of figures, quoted above 
as erroneous, are in agreement. Furthermore, the quotation from Gregory concerning the 
tables is identical with what appears in the first edition. In at least four of the cases cited 
by Mr. Shartle, errors made by Gregory, in seeing one of the tables through the press, are 
definitely suggested. “‘H. Goodwin” is undoubtedly the Henry Goopwyn to whom we 
have already referred, MTAC, v. 1, p. 22-23, 67, 100, 372. According to the Gentleman's 
Magazine, v. 94, p. 381, he died at Blackheath Sept. 27, 1824, aged 79. The most extensive 
tables (in decimal places) listed in the FLETCHER, MILLER & ROSENHEAD, Index, are as 
follows: Circumference, for D = [2(2)200; 6D], in L. Potin, Formules et Tables Numériques, 
Paris, 1925, p. 204-205; Area, for D = [0(.1)100; 6D], in °C. Topp, A Series of Tables. 
Second ed., London, 1853. In Wé6lffing’s Mathematischer Biicherschatz we find °M. F. Kunze, 
7-stellige Kreisflichen von 0.01 bis 99.99. Dresden, 1868; and in Kayser’s Biicher-Lexicon 
is the following: °J. Moprer, Tabellen nebst Anweisung zur bequemen Aufsuchung d. Kreis- 
umfangs, d. Kreisflache, der Oberfliche der Kugel, d. Kubtkinhalis der Kugel bis auf 20 Dezimal- 
stellen, sowie zur Aufsuchg. d. Diameters resp. Radius wenn e. der obigen Dimensionen gegeben 
ist, f. alle Bauhandwerker u. Solche, welche m. Vermessungsarbeiten zu thun haben, sowie zum 
Selbstunterricht . . . Wiesbaden, 1886, 30 p. 

A copy of the first London edition of Gregory’s volume is in the Athenaeum Library, 
Providence, R. I. The first American edition of Gregory (from the second London edition 
of 1833) was published at Philadelphia in 1834. 

In DeMorgan’s article on “‘Table”’ in the English Cyclopaedia, Arts and Science Section, 
v. 7, 1861, we find the following in col. 982: “‘Mr. Goodwyn (of Blackheath) was an inde- 
fatigable calculator, and the preceding tables are the only ones of the kind which are pub- 
lished. His manuscripts, an enormous mass of similar calculations, came into the possession 
of Dr. Olinthus Gregory, and were purchased by the Royal Society at the sale of his 
[Gregory’s] books in 1842.” It would be of interest to learn about the contents of Goodwyn’s 
manuscripts, which are now available at Burlington House. 


82. NYMTP, Table of Characteristic Values of Mathieu's Differential Equa- 
tion, 1945. See MTAC, v. 2, RMT 290. 
For r 2 0, read r > 0, in the following: p. xiii, 1. 10 and 14; p. xiv, I. 5, 10, and 14; 
p. xv, I. 12. 
NYMTP 


83. N. NIELSEN, Recherches Numériques sur certaines Formes Quadratiques. 
Copenhagen, 1931. xi, 160 p. Tables I-II, p. 1-17. 


Errors in T. I-II are indicated in my “A List of errors in a table of numbers D for 


which x? — Dy? = — 1 has solutions in integers,’’ Nieuw Archief v. Wiskunde, s. 2, v. 21, 
1943, p. 194-196. Table I gives all numbers D under 10000 for which the equation 
(1) 2 — Dy 


is solvable; Table II gives those D’s which are sums of two squares and less than 10000 
for which (1) is insolvable. 
The following is a complete errata list. 


Table I. Insert: 2153, 4330, 4777, 5545, 5597, 5669, 6418, 7274, 7330, 8125, 9434, 
9442, 9802. 
Delete: 2157, 2353, 4705, 4786, 5218, 5986, 5989, 6338, 6505, 6605, 7294, 
8585, 8725, 9477, 9962. 
Table II. Insert: 2353, 4546, 4705, 5218, 5986, 5989, 6338, 6505, 6605, 8585, 8725, 
9701, 9962. 
Delete: 4330, 4777, 5545, 5597, 6418, 7330, 8125, 9442, 9802. 


This list was obtained by a recalculation and comparison with two other tables: 


P. SEELING, ‘‘Ueber die Auflésung der Gleichung x? — Ay? = + 1 in ganzen Zahlen, 
wo A positiv und kein vollstandig Quadrat sein muss,’’ Archiv Math. Phys., v. 52, 1871, 
p. 48-49. 
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°W. Patz, Tafel der regelmassig Kettenbriiche fiir Quadratwurzeln aus den natiirlichen 
Zahien 1-10000. Leipzig, 1941. Errata were found in these tables as follows: 

Seeling: Insert 3802; delete 4717. 

Patz: The diamond sign is printed one line too low at 6938, 6949, 6953, 9698. 


N. G. W. H. BEEGER 
Nicolaas Wirsenkade 10, 
Amsterdam 


UNPUBLISHED MATHEMATICAL TABLES 


References have been made to unpublished tables in RMT 279 (Beeger) ; 
MTE 81 (Goodwyn); N 52 (NYMTP); QR 23 (Hillman & Salzer). 


42[F].__N. G. W. H. BEEGER, Second List of Additions to Binomial Fac- 
torisations, Manuscript table of 5 p. in possession of D. H. L. 


The title refers to A. J. C. Cunningham’s extensive Binomial Factorisations in 9 v., 
London, 1923-1929. The first list of additions was prepared by Beeger in 1933 and pub- 
lished under the title Supplement to Binomial Factorisations by the late Lt.-Col. Allan J. C. 
Cunningham, R.E. vol. I to IX, by the London Mathematical Society, London, 1933. The 
present second list gives the factorizations of 81 numbers left unfactored in v. 1 and v. 2 
and the 1933 Supplement. 


D. H. L. 


43[F].—Hansraj Gupta, Three ms. tables in possession of the author at 
Government College, Hoshiarpur, India. 


(a) Tables of partitions giving the values of p(m, m), the number of partitions of m into 
exactly m non-zero summands for values of n < 200. This table gives also the number 
of partitions of m in which the largest part is exactly m. [Compare MTAC, v. 1, p. 
313-314.—Ed.] 

(b) Tables of distributions giving the values of u(m, a), the number of ways in which m 
different particles can be accommodated in exactly a similar cells. m < 50. 

(c) The number of lattice points on the sphere x? + y* + 2? = n, for values of n < 10000. 


H. Gupta 


44[F].—H. K. Hammer, Tables of Periods of Reciprocals of Primes in Various 
Number Systems. Typed mss. prepared by, and in possession of the 
author at 21 West Street, New York City 6. There are copies of these 
mss. in the Library of Brown University. 11 leaves + 35 leaves. 
21 X 29.5 cm. 


Each of the first 11 leaves is devoted to the periods of 1/p, p < 100, in each one of 
the number systems with base 2(1)12. The number of integers in each period, ¢, varies from 
1 to 96. The next 35 leaves contain the periods 1/p, p < 1000 in the decimal notation; 
t varies from 1 to 982. 


H. K. HAMMER 


Epiroriat Note: All of the periods listed on the last 35 leaves were included in those 
published by HENry Goopwyn in his A Table of the Circles arising from the Division of a 
Unit or any other whole number, by all the integers from 1 to 1024, being all the Pure Decimal 
Quotients thai can arise from this source. London, 1823. Tables of the period of 1/p, for 
pb < 1000, were given by K. F. Gauss, in his Werke, v. 2, 1863, second ed., 1876, p. 412-434. 
A table of periods of 1/p to the base 2, for p < 383, has been given by G. BELLAVITIS, in 
R. Accad. Naz. d. Lincei, Cl. d. Sci. fis., mat., e nat., Memorie, s. 3, v. 1, 1877, p. 790-794. 
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A. J. C. Cunnincuam, “On binal fractions,” Math. Gazette, v. 4, 1908, p. 266-267, gives 
for n < 100; and tables to the same extent for the 


bases 3 and 5, in eee Fe etc. fractions,” Math. Gazette, v. 6, 1911, p. 111-116. 
The first 16 digits in the — of 1/n, for m = 1(1)1728, to base 12, may be obtained from 
G.S. Terry, Duodecimal Arithmetic, London, New York, Toronto, Longmans, 1938, p.- 59-70. 


45(G, I].—H. E. Sauzer, Coefficients in Polynomials. Two mss. in the pos- 
session of the author at NYMTP. 


The first ms. gives the coefficients of the various powers of p in the polynomials 
a(qg* — 1*)(g? — 2*) --- (¢ — n*) where g = 1 — p, and m = 1(1)10. 
The second ms. gives the exact values of the coefficients of the eleventh to the twen- 
tieth Laguerre polynomials. 
H. E. SALzer 


Eprtoriat Note: For Laguerre polynomials compare MTAC, v. 2, p. 31. 


MECHANICAL AIDS TO COMPUTATION 


18(Z].—Max Born, R. Firtu & R. W. “A photoelectric trans- 
former,’’ Nature, v. 156, 22 Dec. 1945, p. 756-757. 


The following is the first paragraph of this preliminary account of the transformer, 
from the department of mathematical physics, University of Edinburgh: ‘“‘The ‘Fourier 
transform’ g(y) of a function f(x), usually defined by the integral 


1 
£0) = fle wae, 


plays an important part in many problems of pure and applied physics. It represents, for 
example, the connexion between the intensity distribution of a wave scattered by matter 
of a certain density distribution, which has to be calculated in a number of acoustical 
and optical problems and, above all, in X-ray crystal analysis work. It further allows the 
resolution of a complicated oscillation into a continuous frequency spectrum of harmonic 
oscillations, which is required in many problems of mechanical and electrical engineering. 
It therefore seems of some importance to have an instrument by which the Fourier trans- 
form of a given function can be automatically and quickly produced. We have now suc- 
ceeded in building up an instrument which produced the graph of the function 


= cos (yx + 


on the screen of a cathode ray oscillograph, from a mask cut out of black paper in the 
shape of the graph of the function f(x), or from a record of this function on a plate or 
film in density variation. Obviously two of the functions g’(y) for two values of 6, say, 
& = 0 and 6 = }z, are equivalent to the complex function g(y) when f(x) = 0 for x <a, 


19[Z].—V. Busu & S. H. CALDWELL, “A new type of differential analyzer,” 
Franklin Institute, J., v. 240, Oct. 1945, p. 255-326; 51 figures and 
illustrations. 


The first proposal for a machine to solve differential equations appears to be due to 
Lord Kelvin, in 1876, but the present activity in the construction of such devices was 
initiated by the independent work of V. Bush, and others at the Massachusetts Institute 
of Technology, starting in 1925. The first machine capable of handling a fairly wide class 
of differential equations was described by Bush in 1931, and named the “Differential 
Analyzer.” Since that time many similar machines have been built in this country, in Eng- 
land, Norway and Germany. 
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Engineers and physicists reduce many of their problems to differential equations, and 
families of particular solutions of those equations have proved to be of great service. The 
Massachusetts Institute of Technology and the Rockefeller Foundation joined forces in 
1935 in an effort to produce a more flexible and accurate differential analyzer for the calcu- 
lation of such solutions. The Differential Analyzer, which resulted from this program, was 
placed in service in 1942, and has been in use throughout the war. In November 1945, it 
was opened to the public for scientific work in general including that of industry. 

The natural or canonical form of differential equation for the Analyzer differs some- 
what from the forms usually chosen in the analytic study of differential equations. It may 
be written (the following notation is not used by the authors) 


Ky + 
Xp = + Cp, tr = x + x:+D,, Xu = + Ex, 


where i, j, k, , q, 7, S, t, u, v range over suitable restricted integral values. These equa- 
tions represent the types of mechanism available for the solution. 

The variables, x, are represented in the Analyzer by the angular positions of shafts, 
the shafts being continuously rotated during the solution of a problem. The first equation 
of the canonical set is ‘‘solved’’ by devices called ‘integrators,’ each of which consists of 
a wheel rolling on the surface of a rotating disk. It is interesting to note that the variable 
of integration may be any of the variables in the problem, and in any one problem, integrals 
may be taken with respect to several distinct variables. This fact is frequently relied upon 
to simplify equations. 

The remaining equations of the canonical set are mechanized by means of gear boxes, 
(multiplication by constants Ay¢), “adders,” commonly known as differentials (addition) 
and function tables, on which graphs of the functions f, are plotted. 

The constants, C,, D,, E, and xjo must be set into the machine at the start of each 
problem. Their magnitudes are assigned by the initial conditions of each required solution. 
They are introduced to the machine at scales which are selected with due regard to the 
ranges over which the variables will run during the solution, the ranges available in the 
equipment, and the accuracy with which the solution is desired. 

The M.I.T. Analyzer represents an advance over previous analyzers to some extent in 
its increased accuracy, but to a far greater extent in the speed with which new problems can 
be set up. The accuracy of such complicated mechanisms is difficult to evaluate, even if 
an agreement can be reached on the definition of the term. Errors in any continuous com- 
puter, of which the Analyzer is an example, depend critically on the state of maintenance 
of the parts, on dirt, temperature, wear, etc. etc., as well as upon the wisdom used in the 
choice of scale factors and speeds of solution employed. 

The errors in the solution depend, of course, upon the type of differential equation 
being solved, some equations being sensitive to inaccuracies in the elements and others 
being very stable. 

In view of the facts noted, it is dangerous to quote any figures, but in order to give 
the uninformed a rough idea, it may be said that the errors in a trajectory are likely to be 
of the order of one part in ten thousand, with a very liberal factor one way or the other. 

The speed with which new probiems can be set up on the machine, once they have been 
reduced to the form of codes on punched tapes, is such that only about 15 minutes of 
machine time is required. This contrasts favorably with the one or two days required to 
set up the previous analyzers. 

The class of differential equations which can be handled with the Analyzer is some- 
what wider than that heretofore treated, because the number of computing elements— 
integrators, gear boxes, etc.—is greater. There are 18 integrators available. This does not 
ordinarily mean that 18th order equations can be solved, since some of the integrators are 
likely to be needed for the generation of special functions. For example, if a product xjx 
appears in the original differential equation, it is replaced by 


= f xdat mdx;, 
which requires two integrators to generate. 
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In addition to such functions as may be defined by subsidiary differential equations, 
the Analyzer can accept problems involving three arbitrary functions of one variable each. 
There is no means for dealing with functions of two or more variables, unless these func- 
tions can be reduced to the solutions of subsidiary differential equations, or to functions 
of one variable. There are in existence plans for four “function units’’ which will accept 
values tabulated at discrete intervals and interpolate continuously between them. These 
units are not yet in use. 

Results of the computation are presented either graphically or in the form of typed 
tables. 

Special techniques and schematic symbolisms have been developed for use in reducing 
mathematics to machine connections and are described by Bush and Caldwell. The sym- 
bols are simple and suggest the mechanical units they are intended to represent. While 
not logically necessary, the schematic or short-hand notation is important practically, as 
an aid to the rapid and easy reduction of problems to a form acceptable to the analyzer. 

The design and construction of continuous computers having the accuracy and flexi- 
bility of the M.I.T. Analyzer, demands a high order of ingenuity and mechanical skill. 
Some of the outstanding problems met and solved by the Massachusetts Institute of 
Technology group and their associates are described by Bush and Caldwell. Among these 
are the design of data-transmission and servo units, the development of codes and related 
equipment for the rapid set-up of problems, and the conversion of continuous variables 
(shaft rotations) into digital form, suitable for printing. 

The M.I.T. Analyzer is finding and will continue to find a field of usefulness in the 
solution of the differential equations of engineering and physics, especially as engineers 
delve further into situations which demand the solution of non-linear relations. 


G. R. St1pitz 


20[Z].— INTERNATIONAL BusINESS MACHINES CorP., IBM Automatic Se- 
quence Controlled Calculator, New York, 1.B.M., 590 Madison Avenue, 
1945. 6 p. (21.5 XK 28 cm.) + folding plate (83 X 28 cm.). Copies may 
be procured on application to the publisher. 


This publication contains some notes on the development of this electro-mechanical 
calculator installed and constantly working in the Cruft Memorial Laboratory at Harvard 
University; it was officially turned over to the University on August 7, 1944, and has been, 
and still is, exclusively used by the U. S. Navy. The folding plate gives a complete view 
of the front side of the machine, and various statistics concerning the number of elements 
in the machine and the material used in its construction. 


E. Morre “A slide rule for the addition of squares,” 
Science, n.s., v. 103, 25 Jan. 1946, p. 113-114. 


The author describes a slide rule which he constructed, with appropriate square and 
square-root scales, for evaluating such an expression as d = (x* + ¥* + 2*)}. 


NOTES 
For a Note on Harry Bateman see RMT 289. 


51. Earty DecIMAL DIVISION OF THE SEXAGESIMAL DEGREE (see N 29, 
p. 400f; N 41, p. 454).—In our previous notes on this topic we listed 8 or 
9 editions of De Thiende, 1585, by Stmon STEvIN; but we forgot to give a 
reference to the Dutch edition, 9 or 10, of EzecnreL DE Decker, Eerste 
Deel van de Nievwe Telkonst . . . Noch is hier achter byghevoeght de Thiende 
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92 NOTES 


van Symon Stevin van Brugghe. Gouda, 1626. It occupies the last 27 p. of 
the volume and has its own title page. 


R. C. A. 


52. THE Function LN F (4m, 3, x) AND OTHER RELATED FUNCTIONS.— 
The NYMTP has carried out extensive computations on the function 
In F($n, 3, x) and other related functions to be mentioned below for 
m = 1(2)201 and x = 0(.01).6(.1)2(.2)7(1)45(5)100. The confluent hyper- 
geometric function F(}, }, x) is the solution of the differential equation 


(1) 2xF” + (1 — 2x)F’ — nF = 0, 


satisfying the boundary condition F(}m, 3,0) = 1. If we put F = e!x~tv, 
then it is readily found that » satisfies the differential equation 


16x*0”’ + [3 + (4 — 8n)x — 4x?]v = 0; v(0) = 1. 
For n > x, the Brillouin-Kramers-Wentzel approximation method yields 


v ~ xt-e¥inz, whence F ~ et? ~ ¢Vinz, 


and therefore In F ~ V2nx. The function In F/¥2nx is considerably smoother 
than either F or In F, and for this reason will be tabulated over the major 
range of m and x. 


For x small, neither F nor In F nor In F/V2nx differences well for fixed 


n and running x, but the function F/cosh V(2m — 1)x does; it is therefore 
contemplated to tabulate the latter function for small values of x. That 


F ~ cosh V(2n — 1)x for small x may be shown as follows: If in (1) we 
make the substitution x = y* then 


dF 


If we put F = en, then 

(2) u”+(1—2n—y)u=0; u(0) = 1. 

For x (and therefore y) small, (2) yields 
u’+(1—2n)u=0; u(0) = 1; 


whence u = cosh V2n — 1 y, and therefore F ~ e” cosh V2n — 1y or 


F ~ cosh V2n — 1y. 


NYMTP 


ARNOLD N. Lowan & WILLIAM HORENSTEIN 


53. Guipe (MTAC, no. 7), Suppt. 3 (for Suppl. 1-2, see MTAC, v. 1, 
p. 403-404, v. 2, p. 59).—In MTAC, v. 1, p. 289, is listed C. C. Furnas, 
“Evaluation of the modified Bessel function of the first kind and zeroth 
order,” Amer. Math. Mo., v. 37, 1930, with a table, p. 287, of Jo(x), 
x = [12(2)20(10)100(100)1000; 4S]. On p. 492 of the same volume, S. A. 
SCHELKUNOFF gives a revised form of the table to 5S, from which the Furnas 
values differ by from —0.149 to +0.111 times the corresponding power of 10. 
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Since J. W. WreEncu, Jr. had calculated J)(30) to 20S, MTAC, v. 1, 
p. 200, it was suggested to him that he might check the whole of Schelkunoff’s 
table. He has now (26 Nov. 1945) made the following report: “I recalcu- 
iatted the entire table to 8 and 9S, and thereby proved Schelkunoff’s values 
co be entirely free from error.” 

This calculation also enabled him to check 6 values [x = 12(2)20, 30] of 
I(x) in the 6S table of H. E. H. Wrincn & D. M. Wrincu, Phil. Mag., 
s. 6, v. 45, 1923, p. 847, and v. 47, 1924, p. 63. Compare MTAC, v. 1, p. 226. 
The values for x = 12 and 14 checked exactly; for x = 16 the last digit 
needs to be increased by unity; for x = 18 an increase of 3 units is neces- 
sary; for x = 20, a decrease of 3 units; and for x = 30 a decrease of 2 units. 

R. C. A. 


54. A Huce NumBeEr.—Introductory to comment on a recent paper 
involving tables, I shall reprint some of the items published in my note on 
“Huge Numbers,” Amer. Math. Mo., v. 28, 1921, p. 393-394. What is the 
largest number that we can express by three digits? One answer is 9 or 
g3s7420489, C. A. LAISANT drew attention to this number in his Initiations 
Mathématiques, Paris, 1906 (English edition, London, 1913, p. 101). He 
then remarks that in decimal numeration this number would have 369693100 
figures. To write it on a single strip of paper, supposing that each figure 
occupied a space of one fifth of an inch, the length of the strip would be 
1166 miles, 1690 yards, 1 foot, 8 inches. 

Writing in 1913 A. C. D. CROMMELIN stated, Br. Astr. Assoc., J., 
v. 23, p. 380-381, that he had come across the problem “‘in an old log- 
arithm book.” By the aid of 61-figure logarithms Crommelin found 
log N = 369,693,099.6315703587 ---; whence the number of figures indi- 
cated above. He then found the first 28 figures of N to be 428,124,773, 
175,747 ,048,036,987,115,9, and the last three to be 289. Crommelin then 
continues, ““A knowledge of 30 figures out of 300 million may seem trifling, 
but in reality the error involved in taking all the remaining figures as zeros 
is only one part in a thousand quadrillions.' If the number were printed 
with 16 figures to an inch (about the lightest packing for decent legibility), 
it would extend over 364.7 miles. . . . If printed in a series of large vol- 
umes we might get 14,000 figures to a page, and with 800 pages to the volume 
it would fill 33 volumes. There are more than twice as many digits in the 
number as there are letters in the whole of the Encyclopedia Britannica. .. . 
We may take it as morally certain that we can write with three digits a 
number vastly exceeding the number of electrons in the whole of creation, 
which is a somewhat startling fact. Indeed, even the number 4“ (which is 
13407813 followed by 147 other figures) probably exceeds the number of 
electrons in creation.” 

In Br. Astr. Assoc., J., v. 26, 1915, p. 46-47, D. G. McINTYRE gave 
the last eight figures as 17,177,289. 

We now turn to Matematisk Tidsskrifi, A, 1941, p. 63-70, and an article 
by Cur. WEIss, entitled “ ‘Hu’, Tallet 9 og Endecifrene i Potenser 
af 9.” Weiss found log N to 74D as follows: 369693099.631570 358743 
543095 095482 929683 400246 855547 962405 328896 342699 975257 003407 
00; whence he derived the first 60 figures of N to be (details of computation 
given) 4281 247731 757470 480369 871159 305635 213390 554822 414435 
141747 53. 
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94 QUERIES—REPLIES 
i On p. 67 is a table of the last 26 figures for each of 35 values of 9, logai 
E n = 89,100,200, - - -,387420489. Thus the last 26 figures of N are found intes 
to be 
: 24 178799 359681 422627 177289. Ms 
2 These results check with those quoted above, except in the case of the _ 
* first of the McIntyre figures. Weiss gives also two tables and formulae for us. F 
wv finding last figures of 9". only 
; J. W. MearEs in Br. Astron. Assoc., J., v. 31, 1921, p. 277-278, com- 
ments on 9!" and finds that its value is greater than 10 to the power 2 
107999 but less than 10 to the power 107001, = 
R.C.A. 
heat 
1In accordance with British usage, Crommelin here means 1000 X 10%; in the United I hz 
States this would be interpreted as 1000 x 10%. argt 
m 
55. A NEw REsuLT CONCERNING A MERSENNE NUMBER.—(Compare ‘a 
N. 23, 33, v. 1, p. 333, 404). On 9 February 1946 I finished testing the Br 
character of the Mersenne number Moe» = 2”° — 1 = 8627 18293 34882 Re 
04734 29344 48278 46281 81556 38862 15212 98319 39531 55279 74911. 
Since the final residue, the 228th, was not zero the conclusion is that Ma, {| y 
is composite. 336: 
The Lucasian sequence used was 4,14,194,37634, 1416317954, etc. Equ 
: The 228th residue was found to be 1970 11660 94225 75309 56180 (see 
: 91126 86257 27776 96596 41856 06805 84362 68648 91891. root 
: Thus, among these numbers M,, up to and including p = 257, there tan 
are only three whose characters are unknown, namely: p = 193,199,227. tion 
There are, however, eleven M,, known to be composite, but of which no of c 
factor is known. fron 
I have begun a similar investigation of M19. the 
Horace S. UHLER oth 
12 Hawthorne Avenue inve 
Hamden 14, Connecticut is a 
QUERIES 
17. TABLES FOR CircLEs.—In O. G. GreGory, Mathematics for Prac- 
tical Men, London, 1825, p. 406, after “A Table of Circles, from which 
knowing the diameters, the areas, circumferences, and sides of equal squares P. 2 
are found,”’ by Goopwyn (see MTE 81), Gregory remarks that this table P. 2 
was ‘‘to supersede the necessity of consulting some erroneous tables of the P.2 
areas, &c. of circles recently put into circulation.” What author and publi- P.2 
cation are here indicated? The English Catalogue lists the following anony- P.2 
mous item issued in the following year: Tables of Areas and Circumferences P.2 
of Circles, 3 parts, London, 1826. P.2 
a. P.2 
QUERIES—REPLIES P. 2 
21. Briccs’ ARITHMETICA LOGARITHMICA (Q7, v. 1, p. 170).—In my P. 2 
library is a copy of this volume with the extra 12 pages containing the P. 2 
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logarithms of numbers [100001(1)101000; 14D] and the square roots of 
integers [1(1)200; 11D], with first differences in each case. 
ALBERT E. LOWNES 


16 Barberry Hill 

Providence, Rhode Island 
EpitoriaL Note: At present only four existing copies of these extra are known to 
us. Rather extensive inquiry suggests that the copy owned by Mr. may be the 
only one in America. 


22. INTEGRAL AND FUNCTIONAL TABLES.—(QI15, p. 459).—With refer- 
ence to tables of foe-“dt, I should like to draw attention to my paper 
“An approximate method for calculating heat flow in an infinite medium 
heated by a cylinder,’’ Phys. So. London, Proc., v. 56, 1944, p. 365, where 
I have given a rough chart of this integral [24f,""e-“dt] for a complex 
argument [8 = 0(15°)90°; w = 0(.1).3(.2).7, 1], together with formulae for 
computation in a number of cases. 

S. WHITEHEAD 
British Electrical and Allied Industries 
Research Assoc., London, England 


23. Roots OF THE EQUATION TAN x = cx (Q8, v. 1, p. 203; QR10, p. 
336; 17, p. 459).—The NDRC report, Tables for Solutions of the Wave 
Equation for Rectangular and Circular Boundaries having Finite Impedance 
(see MTAC, v. 1, p. 438-440) has (in a different notation) the first four 
roots (i.e. the first four branches of the functions) of the equations 
tan x = cx, tanh x = cx, cot x = cx, and cothx = cx, tabulated as func- 
tions of the complex variable c. We have obtained expansions of x in powers 
of c and 1/c (used in preparation of the above mentioned NDRC report) 
from which one can easily compute values on the other branches or improve 
the accuracy of the above tables. We also have theoretical material and 
other types of expansions for the regions about the singularities of the 
inverse functions. This material is in our possession at the NYMTP, but 
is as yet available only to certain Government agencies and activities. 

A. P. Hittman & H. E. SALZER 
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P. 205, a 17, for je, read 
P. 213, As 38, 1. 3, add x = 0(.5)12; Bz 4, for 5D, read 4D. 
P. 214, B; 13, for read and for 4, read 4; B; 17, for read 
P. 215, B; 8, add also 5D, = 1(1)8; Bz 22, for =, read or. 
. 217, Ce. 7, delete = 14D; D, 7, 1. 3, when m = 10, 11, read s = 3, 3. 
221, As 1, delete x } 21.6; As 1, for x } 21.6, read for each n; As 6, 7, for Yo(x), Yi(x), 

read Y(x), Y™(x), and transfer these two entries as C; 8, 9. 
222, As 12, for zeros, read first zero. 
226, Ai 12, 13, for 5D, read 5S, and for 4D, read 4S. 
227, As 28, for OKayYaA 2, x, read OKAYA 2, x71 
241, B 5, for read J_1;6(4x+); and for 41.035, read 41.305. 
245, 1. 14, for E,(x), read E,(x). 
251, B 16, read 4D, Ki2(x), SEELIGER, x = .25(.01)1(.05)7. 
284, 1. 13, for Yn(x), for, read Y,(x) for each n, for. 


nd 
the 

for 
m- 

ver 

ited 

are 

the 

382 

11. 

ere 
27. 
no 

4 

ac- 

ich 
res 
ble 

the 

ces 


P. 


. 6, 1. —8, for read 


96 
P. 


289, 1. 4, for 393, read 395. 
414, 1. 8, for 1942, read 1941. 
476, Riche de Prony, for Caspard, read Gaspard. 


V.2 


1, L —21, for St., read 51. 
5, 1. 1, for Tablitsy, read Tabliisy. 


6 
. 9, 1. 11, for z or, read or z,. 
1 
1 
1 


6, 1. 32, 4.7, for Sums of Products, read Sums and Products; 1. —11, for Gennochi’s, 
read Genocchi’s. 


29, 1. 32, for give complex roots, read give complex roots directly as decimals. 
60, 1. —1, for arccos(sinhy), read arccos(—y/sinhy). 
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